Simulation of induction at low magnetic Prandtl number 
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We consider the induction of magnetic field in flows of electrically conducting fluid at low magnetic 
Prandtl number and large kinetic Reynolds number. Using the separation between the magnetic 
and kinetic diffusive lengthscales, we propose a new numerical approach. The coupled magnetic 
and fluid equations are solved using a mixed scheme, where the magnetic field fluctuations are fully 
resolved and the velocity fluctuations at small scale are modelled using a Large Eddy Simulation 
(LES) scheme. We study the response of a forced Taylor-Green flow to an externally applied field: 
tology of the mean induction and time fluctuations at fixed locations. The results are in remarkable 
agreement with existing experimental data; a global 1/ f behavior at long times is also evidenced. 

PACS numbers: 47.27.Eq,47.65+a,52.65Kj,91.25Cw 



One of the strongest motivation in the study of non- 
linear effects in magnetohydrodynamics is that electri- 
cally conductive flows are capable of dynamo action: the 
stretching of magnetic field lines by the flow velocity gra- 
dients can exceed the (Joule) diffusion. A bifurcation 
threshold occurs, above which the self-generation of a 
magnetic field takes place. It has been validated in con- 
strained flows of liquid sodium which mimic analytical 
models: the Karlsruhe [l| and Riga experiments p]. The 
self-generation of a magnetic field in non-constrained ho- 
mogeneous flows is still an open problem actively studied 
by many groups In this research, numerical stud- 
ies have long played an important role. Kinematic dy- 
namo simulations assume a given pattern of a station- 
ary velocity field and study the initial linear growth rate 
of magnetic field perturbations. They have been exten- 
sively used to test the dynamo capacity of fiow geometries 
and proved to be successful at determining the dynamo 
threshold in the Karlsruhe and Riga experiments 0, 0| . 
They have also shown that dynamo action is a possibility 
in unconstrained homogeneous flows of the von Karman 
type 1^ 0|- Another numerical approach is to perform 
Direct Numerical Simulations (DNS) of the full govern- 
ing equations: the induction equation coupled with the 
fluid dynamical one by the Lorentz force, the flow be- 
ing sustained by a given force (or equivalently an aver- 
age geometry). They have confirmed that dynamo ac- 
tion is present in flows with differential rotations and 
helicity 0, 0, 0|. However, DNS are at present re- 
stricted to situations where the magnetic Prandtl num- 
ber, P„i = v/X (where A is the magnetic diffusivity) is 
of order one, i.e. to situations where the smallest scales 
of the magnetic and velocity fields have the same char- 
acteristic size This is not the case in liquid metals, 
which have very small magnetic Prandtl number values 
e.g. Pm ~ 10~^ for hquid Gallium and Pm 10~^ for 
liquid Sodium. Recall that, below the dynamo threshold. 



a stationary forced fiow with a power input e (in Watts 
per kg) has a viscous dissipative scale ry„ ^ {v^ jt)^^'^ and 
a Joule diffusive scale r\B ~ (A'^/e)^/^ — hence a ratio 
tIuMb ^ Pnl^ ■ Therefore, at low P^, the magnetic dif- 
fusive length scale is very much larger than the velocity 
dissipative scale. If non-linear effects are to develop, the 
magnetic Reynolds number Rm ~ UL/\ (where U and 
L represent the characteristic velocity and scale of the 
flow) must be at least of order one and thus the kinetic 
Reynolds number of the flow. Re ~ UL/v ~ Rm/Pm, 
must be very large (turbulence is fully developed). A 
DNS treatment of such a system is at present out of reach. 

In this paper, we present a new approach for the study 
of the magnetic induction in large Re - low Pm flows; we 
restrict ourselves to regimes below the dynamo threshold. 
In this parameter region, the magnetic field "lives" essen- 
tially within the large and inertial hydrodynamic scales. 
We thus propose to treat with a sub-grid model the veloc- 
ity scales which are smaller than the magnetic diffusive 
length. Schemes using hyperviscosity have previously 
been used Here, we prefer a LES approach, which 

has proved very successful for the simulation of turbulent 
fiows with large scale structures and for the modelling of 
energy transfers this hybrid scheme, we solve 

the induction equation on a fully resolved grid and we 
use a LES method for the velocity field, with a cut-off 
scale at the end of the magnetic diffusive range. We con- 
sider the response of a conductive fiuid to an uniform 
magnetic field: topology of the mean induced field and 
spatio-temporal features of the magnetic fluctuations are 
studied. The chosen flow is a forced Taylor-Green vortex 
(TG). It shares many similarities with the experimen- 
tal von Karman swirling flows which have already been 
investigated in DNS near P,„ - 0(1) ji,E3- 

In non-dimensional form, the incompressible MHD 
equations have two classical control parameters, the mag- 
netic and kinetic Reynolds numbers, and one has to 
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choose a forcing mechanism that generates the desired 
values of Rm and Re. In order to be closer to experi- 
mental procedures, we prefer to fix the driving force and 
the magnetic Prandtl number. Hence, the dynamical 
time to is set to the magnetic diffusion time scale, i.e. 
toX/L'^ ~ 0{1), where L is a length scale characteristic 
of the system size. Changes in magnetic diffusivity for 
real fluids would change that time scale. We write the 
MHD equations, with constant unit density, as 



dtu + u.Vu 
dth 
V.u 



V X (u X B) + V^b , 
, V.b = , 



(V X b) X ^1) 

(2) 
(3) 



where u is the velocity field, B = Bq + b is the net 
magnetic field in the flow, sum of the applied and in- 
duced fields. Once the amplitude F of the driving force 
is fixed, the (non-dimensional) rms intensity of the ve- 
locity fluctuations is Urms ~ the Reynolds num- 
ber is Re \/F I Pra and the magnetic Reynolds num- 
ber is Rm ^ ypF . When the interaction parameter, ra- 
tio of the Lorentz force to the inertial forces, defined as 
iV ~ Bf^jurms ~ Bq/^/F is small, the back reaction 
of the induced field on the velocity field is negligeable. 
The above expressions are only dimensional estimates; in 
practice, the characteristic flow quantities are computed 
as mean temporal values from the data - cf. Table 1 . 

We use a parallelized pseudo-spectral code in a [0— 27r]'^ 
periodic box. Time stepping is done with an exponential 
forward Euler-Adams-Bashford scheme. The LES model 
is of the ChoUet-Lesieur type in which the kinematic 
viscosity u is replaced in spectral space by an eddy vis- 
cosity. In Eq. the magnetic Prandtl number is then 
replaced by: 

Pm{k,t) = Q.l{l + b{k/K,f)y/E,{k^K,,t)/K, . (4) 

Here is the cut-off wavenumber of the velocity field, 
and Ev{k,t) is the one-dimensional kinetic energy spec- 
trum. The effective Prandtl number Pm^a is obtained as 
the temporal mean of Pm{0,t). Note that the effective 
fluid viscosity v^s is of the same magnitude. A consis- 
tency condition for our approach is that the magnetic 
field fluctuations are fully resolved when 2tt/Kc is smaller 
than the magnetic diffusive scale jys ~ lo/Rm^, Iq being 
the integral scale computed from the kinetic energy spec- 
trum. The flow is driven by the TG vortex geometry 



' TG 



(fco) = 2F 



sin(fco x) cos(fco y) cos(fco z) 
— cos(fco x) sin(fco y) cos(fco z) 




(5) 



(fco, ko, fco) is the wavevector that prescribes the velocity 
large scale (hereafter ko = 1). The Ftg and Bq am- 
plitudes are chosen such that the interaction parameter 
N remains smaller than 10~^. After an initial transient 
{t < 10) the flow has reached a steady state: the kinetic 



energy fluctuates less than 3.5% around its mean value. 
All quantities are tracked up to ^max = 410^0 — note 
that 200to is of the order of the measurement time in 
most Sodium experiments |l4l llfj . Ilfij ] . For comparison, 
the eddy turnover time ttvl ~ lo/urms is given in Table 
1. 
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feo = 1 


Rm = 6.65 


Rm = 6.68 




Hij, ~ yo.94 


Hij, = yo.yo 


128^ grid points 


P™„,,. ~ 7.22 10-* 


Pm^a = 7.26 10"" 


— ^max 3 


N = 8.23 10-^ 


TV = 8.18 10"^ 


fcmax — 64 


^0 = 2.338 


lo = 2.337 


tmax = 410 


It = 0.024 


It = 0.024 




riB = 0.565 


77s = 0.563 




TNL = 1.217 


TNL = 1.224 




Urms = 2.843 


'^rTYis — 2.858 




brms — 0.061 


brms = 0.064 




max|u| = 8.211 


max|u| = 8.249 




max|b| = 0.160 


max|b| = 0.180 



TABLE I: Time averaged quantities: 

,2\l/2 



(U^) , br 



b ) , flow integral scale Iq = 27r _E„(fc)/fc/ i5„(fc), 

— 1/2 

loRe , diffusive scale r/B and eddy 
Non-dimensional parameters: effec- 



Taylor microscale It 
turnover time tnl ■ 



tive Prandtl number Pm^ff, kinetic Reynolds number Re 



loUrma / i^cft (sBC text), and magnetic Reynolds number R 
Pm^ffRe, Taylo-based Reynolds number Ri^, 
tion parameter A'^ = RmBo/urms- 



rI^^ , interac- 



Figure 1 shows the power spectra of the velocity and 
magnetic field fluctuations with Bq applied along the 
X-axis (a direction perpendicular to the rotation axis 
of the counter-rotating eddies of the TG cells). The 
kinetic energy spectrum exhibits the k^^/^ Kolmogorov 
scaling law maintained throughout the range by the 
LES scheme. The peak at low wavenumber is due the 
large scale TG forcing, also visible on the magnetic field 
spectrum. The magnetic inertial range is well fitted by 
a fc-^^/"^ power law in agreement with a Kolmogorov 
phenomenology 0, . The magnetic diffusive scale is 
reached within the computational box. The main goal 
of our numerical strategy is thus achieved: the magnetic 
fluctuations are fully resolved in a range of scales 
at which the velocity field follows the Kolomogorov 
self-similar structure of turbulence. Hence, we get 
the possibility to study magnetic induction in a fully 
developped turbulent flow at low magnetic Prandtl 
number. 

Figure 2 displays isosurfaces of the local induced mag- 
netic energy {Eb(yi,t))rp averaged in the time interval 
T = [10 — 410], shown at 80% of its maximum value. For 
comparison, we also plot isosurfaces of the induced mag- 
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FIG. 1: Magnetic (solid line) and kinetic (dash line) energy 
spectra computed at t = 210 for RUN 1 with Bo = 0.1 x. 



nctic energy, {Eb^iini'^,t))rp, obtained numerically from a 
linear approximation based on time averaged velocities: 
AV^b = -Bo V(v(x,t))y. This is similar to numeri- 
cal studies based on the averaged flow geometries 0, ^3 ■ 
When Bq is applied along z, in a direction parallel to 
the rotation axis of the TG eddies, the most intense 
magnetic energy structures are concentrated round the 
z = 7r/2,37r/2 planes, in agreement with the differential 
rotation of the TG vortex. Moreover, the most intense 



structures of (i?;,(x, t))^ and {Eb^un{^,t))j, fields coin- 
cide. For Bo along the x-axis, one observes the main 
induction concentration around the z = 0, tt planes, as 
expected from a direct inspection from the flow forcing. 
However, the most intense structures of the (£^t,(x, i)) 
and {E},^iin{yi,t))j, fields do not coincide everywhere in 
that case (see location (7r/2, 7r/2, 0) in Fig. 2(bottoni), 
for example). Note also that the linear calculation over- 
estimates the time averaged magnetic fluctuations, what- 
ever the orientation of the applied field. Altogether it 
shows than one should be cautious when using average 
velocity fields in the calculation of magnetic induction, 
particularly if restricted to linear effects. The difference 
between the fields is probably linked to the large scale 
electromotive force due to turbulent motions. The in- 
fluence of this force, as well as the large scale induction 
topology and its connection with the small scale fluctua- 
tions, will be reported in a forthcoming paper [20l |. 

Figure 3 shows the temporal fluctuations of the in- 
duced fleld amplitude, |b(x, i)|, probed inside the flow at 
two locations chosen from the previous topological obser- 
vations, for Bo along the x-axis. This is equivalent to us- 
ing local probes as in laboratory experiments. The inten- 
sity of the induced magnetic fleld has strong local fluctu- 
ations. The point at (0, tt, 0) is in a region of strong mean 
induction, whereas the point at (O.Gtt, O.Gtt, O.Gvr) is at lo- 
cation of low mean induction (cf. Fig. 2(bottom)). We 
observe that, occasionally, the induced fleld gets larger 
than the applied fleld. In fact, if small amplitude fluc- 
tuations (about 10%) are induced over time intervals of 
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FIG. 2: Topology of the local induced magnetic energy, av- 
eraged in time, when Bo is applied along the z-axis (top) 
and along the x-axis (bottom) - in red: (i?i,(x, t)) j,; in blue: 
{Eb,im{Ti,t))j, - (see text). The isosurfaces are plotted at 80% 
of the maximum values of the fields : ma,x{Eb)j, = 0.0056 
and ma,x{Eb,un}j' ~ 0.0063 for Bo = 0.1 z, and ma.x{Eb)j, = 
0.0041 and max(£b,Hn)T = 0.0063 for Bo = 0.1 x. 



the order of the diffusive time to, much larger variations 
(~ 300%) can be observed over long time periods, of the 
order of lOio- These observations are in excellent qual- 




FIG. 3: Time traces of |b(x,t)|, for Bo = 0.1 x, at two 
fixed points. In blue: (0,7r,0), mean value {\h{x,t)\}T/ Bo = 
0.92, fluctuation level \h{x,t)\rms/ Bo = 0.28. In red: 
(0.67r, 0.67r, 0.67r) mean value (|b(x, t)|>T/Bo = 0.44, fluctua- 
tion level |b(x, t)|rms/Bo = 0.19. 
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itative agreement with the experimental observations at 
comparable Rm and Pm 0, E E, E| ■ In order to be 
more quantitative, we analyze the time spectra; we fo- 
cus on the case with Bg applied along the x-axis, but 
the results are identical when Bq is along z. We plot 
in Figure 4 the power spectra of the temporal fluctua- 
tions of the magnetic field component (x, t) recorded 
at (0,7r, 0). The higher end of the time spectrum follows 
a behavior close to /~^^^'^, as can be expected from the 
spatial spectrum using the Taylor hypothesis of "frozen" 
field lines advected by the mean flow 0. In addition, 
for frequencies roughly between 1/to and l/lO^o, the time 
spectrum develops a 1 / f behavior, as observed in experi- 
mental measurements . It is not present on the spatial 
spectrum in Figure 1, and thus appears as a distinctive 
feature of the time dynamics of the induced field. It is 
also independant of dynamo action, as it is also observed 
in the Karlsruhe experiments . Finally, our numerical 
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FIG. 4: Power spectral density of the magnetic field fluctua- 
tions of bx{'>i,t) in time, recorded at space location (0,7r,0), 
when Bo = 0.1 x. (a) PSD computed as averages over Fourier 
transforms calculated over long time intervals (~ 164fo) to 
emphasize the low frequency behavior; (b) PSD estimated 
from Fourier transforms over shorter time intervals (~ lOto). 
The behavior is identical for the by (x, t) and bx (x, t) field 
components. 

study reveals one remarkable feature: the 1// behavior is 
a global feature. It is observed on the fluctuations of the 
magnetic energy, as shown in Figure 5 (as a scaling 
regime). We thus propose that it results from induction 
processes which have contributions up to the largest scale 
in the system. 
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FIG. 5: Power spectral density of the time fluctuations of the 
magnetic energy Eb{t) — (b^(t))/2, intergrated over space. 

To summarize, the mixed numerical scheme proposed 
here proves to be a valuable tool for the study of mag- 
netohydrodynamics at low magnetic Prandtl numbers. 



We have considered here the response to an externally 
applied field. The time behavior of magnetic field fluctu- 
ations is found in excellent agreement with experimental 
measurements. It has also revealed that the 1/f regime 
detected locally traces back to the global dynamics of 
the flow. Future work will analyze the contribution of 
turbulent fluctuations to the large scale magnetic field 
dynamics, and the influence of the magnetic Prandtl 
number on the threshold of the dynamo instability. 
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